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We present a detailed analysis of large scale simulations of avalanches in the 2D Abelian sandpile 
model. We compare statistical properties of two different decompositions of avalanches into clusters 
of topplings and waves of topplings. Auxiliary critical exponents are introduced and the existence 
of the exponent governing the contraction of avalanches claimed in our previous work [Priezzhev et 
al, PRL 76,2093 (1996)] is confirmed. We also give more elaborated argumentation for the exact 
values of the exponents characterizing the statistics of waves. 
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I. INTRODUCTION 



The sandpile model introduced by Bak, Tang and 
Wiesenfeld jlj] serves not only as a lapidary formula- 
tion of basic principles of self-organized criticality (SOC) 
but also seems to be an appropriate candidate for exact 
determination of all important critical exponents. In- 
deed, the first steps following Dhar's discovery of the 
Abelian structure of the sandpile model were encourag- 
ing. They include determination of the total number of 
allowed configurations in the recurrent state Q], evalu- 
ation of the height probabilities J|,|| and height-height 
correlation functions PP], interpretation of the inverse 
Laplacian operator A as an expected number of top- 
plings at a given site due to a particle added to another 
one g]. Nevertheless, all analytical results obtained up to 
now catch either static properties of the recurrent state 
or diffusion-like dynamics of individual particles. The 
avalanche dynamics as such, responsible for SOC, slips 
off the analytical description even in the simplified lim- 
iting case of large avalanches. The main obstruction is 
that existing renormalization group methods |6| neglect 
essential peculiarities of the toppling process, and the 
complicated spatio-temporal structure of avalanches pre- 
vents exact evaluation of the critical exponents. 

To advance the analysis of avalanche dynamics, vari- 
ous decompositions of avalanches in the Abelian sand- 
pile model (ASM) into more elementary objects have 
been proposed. In particular, Grassberger and Manna 
noticed that each avalanche can be represented as a 
set of embedded clusters of sites related to a given num- 
ber of topplings. To make use of this construction for 
determination of critical exponents it is desirable to ob- 



tain a dynamical procedure which naturally divides the 
avalanche into a collection of clusters. It means that due 
to the Abelian property of toppling operators, one can try 
to change the order of topplings so that each avalanche 
would expand to the largest cluster and then contract by 
smaller and smaller sets of toppling sites. In our previ- 
ous works g-[l0| we made such an attempt proposing a 
decomposition of avalanches into waves of topplings. 

The main feature of the wave structure of avalanches 
is a possibility to set up a one-to-one correspondence be- 
tween waves and two- rooted spanning trees ||. Using 
the spanning tree representation for waves, one can ap- 
ply the methods of graph theory to calculate the critical 
exponents of wave and avalanche statistics (|Jl^l . 

On the other hand, it has been found out |IT|]§|| that 
the set of waves and the set of clusters for a particular 
avalanche do not coincide. Namely, waves have such an 
irrelevancy in their superimposing that the next wave can 
overlap the previous one and the package of waves does 
not form embedded sets of sites like clusters. The obser- 
vations of Dhar and Manna JTT[ and our simulations on 
small lattices raise hope that the overlapping^ of waves 
are relatively rare events. We have assumed [|10| that one 
can neglect the difference between clusters and waves of 
topplings and consider each next wave embedded into 
the previous one as a typical situation. Based on this 
assumption we suggested a method of evaluation of the 
basic critical exponents of 2D ASM. Our latest simula- 
tions, however, have shown that the next wave typically 
overlaps the previous one. Moreover, the large-scale sim- 
ulations of Paczuski and Boettcher state that the 
average difference of sizes between two subsequent waves 
is actually negative. 
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Nevertheless, we will show in this paper that it is pos- 
sible to modify our simplified scenario of the avalanche 
process and to describe the phases of expansion and con- 
traction in terms of the wave decomposition. Besides, we 
will demonstrate that the theoretical predictions are in 
complete agreement with the numerical data obtained by 
Paczuski and Boettcher |Q . 

The paper is organized as follows. In Section II we 
formulate the ASM model, define avalanche clusters and 
waves of topplings and introduce the basic ideas of ex- 
pansion and contraction of avalanches. Section III is de- 
voted to the description of the local dynamics of waves 
and proof of the existence of the contraction exponent. 
In Section IV we present analytical derivation of the ex- 
ponents of conditional distribution of waves obtained by 
Paczuski and Boettcher |^2| from extensive numerical 
simulations. In Section V, an elucidative point of view 
on the renormalization group approach || to the sandpile 
model is suggested. 



II. BASIC CONCEPTS 

We remind the definitions of the model, waves of top- 
plings, avalanche clusters and explain our basic ideas. 

In 2D ASM one starts from the empty square lattice 
(occupation numbers Zi — for all sites) and drops sand, 
particle by particle, at random sites: Zi — > Zi + 1. If any 
Zi > 4, the site i is unstable and topples: Zj — > zj — Ay 
where A is the Laplacian matrix. The toppling at i may 
cause instability at its nearest neighbors. The subsequent 
topplings continue until there are no more unstable sites. 
Then, one adds again a particle at random site initiat- 
ing a new chain of topplings and so on. The process of 
toppling during each perturbation is called an avalanche 
and the set of toppled sites form a compact cluster of all 
toppled sites. 

To obtain the wave decomposition of an avalanche || , 
one has to topple all sites that become unstable after 
adding a particle at i keeping this site out of the second 
toppling. The set W\ of toppled sites is the first wave 
of topplings. All sites except maybe the site i become 
stable after the first wave. If the resulting height Zi > 4, 
one topples the site i the second time and continues the 
relaxation process not permitting this site to topple the 
third time. The new set W2 of relaxed sites is the sec- 
ond wave. The process continues until the site i becomes 
stable and the avalanche stops. 

Grassberger and Manna defined clusters C n ,n = 
1, M of sites toppled not less than n times during the 
given avalanche, M is the number of topplings at the ini- 
tially perturbed site. The sets C n ,n < M are all compact 
and each C n contains the clusters C ra +i, Cm- 

It is possible to evaluate the asymptotics of cluster 
size distribution considering the set of generated clusters 
without reference to a particular avalanche they belong 
to. According to [0, the expected number of topplings 



at site j due to adding a particle at site i is given by 
the lattice Green function Gij = [A - The number 
of topplings Gij coincides with the expected number of 
clusters C n containing the site j in an avalanche started 
at i. Therefore, the probability that the linear extent r 
of a cluster exceeds the distance \i — j\ between i and j 



Prob(r > \i-j\) ~ G; 



(1) 



Using the known asymptotics of the Green function for 
large distances G(r) ~ ln(r) and compactness of clusters 
(cluster area s c ~ r 2 ), we get 



P(s c ) = P(r) 



dr 

dSr. 
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(2) 



It was established in S, that every wave is a com- 
pact set without holes and each site in a wave topples 
exactly once in that wave. Thus, the expected number 
of topplings Gij given by the lattice Green function can 
be expressed alternatively by the probability that a wave 
taken from an arbitrary avalanche initiated at site i cov- 
ers site j. Writing Eq.( |l|) for waves, we get again the 
size distribution similar to Eq.(0) 



P{s w ) ~ —, 

Sin 



(3) 



where s w is the area of a wave. 

To find critical exponents characterizing the size distri- 
bution of avalanches, one also needs a general picture of 
the avalanche process as a whole. In the case of clusters, 
the picture is quite clear. The set of clusters is ordered 
and each next cluster is embedded into the previous one. 
However, the clusters of topplings, being convenient for a 
computer decomposition of avalanches, are hardly repro- 
ducible by dynamical rules as each cluster grows mono- 
tonically during the whole avalanche process. 

On the contrary, the wave construction admits a sim- 
ple dynamical interpretation but loses the property of 
ordering which is inherent in the case of clusters. In 
spite of the irregularity of waves, we are still able to 
use a partial ordering of waves assuming that a typi- 
cal avalanche consists of two phases: fast expansion and 
slow contraction. The first phase contains relatively few 
waves with a large negative difference between subse- 
quent waves A(sfe) = Sfe — -Sk+i- The second phase forms 
the main body of an avalanche with a positive average 
difference A(sfe) > 0. In [[h| the fast phase was asso- 
ciated with the single first wave which reaches at once 
a maximal size the given avalanche spreads. The posi- 
tive difference for the rest of the waves was assumed to 
be dependent only on the size of a preceding wave and 
satisfied the scaling law 



(AGO) 



(4) 



for large s. 
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FIG. 1. Sizes Sfe of waves in a typical avalanche on the lat- 
tice of size L 2 , L — 500 (empty squares) ; the same quantities 
subtracted by the size of overlapping: Sk — A~(sk-i) (filled 
diamonds) . 

If the law Eq.(|J) is valid for clusters of topplings as 
well, the density of clusters can be defined as the average 
number of clusters of size between s c and s c + ds c in one 
avalanche of the size S > s r 



dn 
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(5) 



By assumption, the density depends on s c but not on 
5. Then, the critical exponent r in the distribution of 
a number of sites covered by an avalanche P(S) ~ S~ T 
can be related with the exponent a in Eq.(0). Indeed, 
the probability distribution of cluster sizes P[s c ) is pro- 
portional to the probability of avalanches whose size S 



exceeds s c : P(S 
clusters Eq.(||): 



> s, 



■) 



-T+l. 



and to the density of 



Comparing Eq.(^) with Eq.(||) we obtain 
r + a = 2, 



(0) 



(7) 



and the problem of finding the basic exponent t is re- 
duced to search for the exponent a which is related more 
directly to details of the avalanche process. 

The 'contraction' exponent a is well defined for 
avalanche clusters or for waves provided that one can ne- 
glect the differences between these two kinds of objects. 
In this connection, the following questions arise. Is it 
possible to define the 'contraction' exponent for waves 
taking into account overlappings and if so, what is a cor- 
respondence between its numerical values for clusters and 
waves? Can we establish the same relation Eq.(^) for 
waves and use their spanning tree structure to estimate 
the critical exponents of the model? 

In the following sections we discuss these questions us- 
ing large scale simulations of clusters of topplings, a more 
elaborated analysis of waves and new numerical data for 
subsequent waves obtained by Paczuski and Boettcher 
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FIG. 2. The average values of A~(s), A + (s) for waves and 
A(s) for clusters as functions of their size s (see text for defini- 
tions), obtained from the simulations data of 10 6 avalnches on 
lattice of size L — 500. The graph A(s) is shifted vertically. 



III. CONTRACTION EXPONENT FOR WAVES 



First, we present a picture of a typical avalanche of 
2D ASM on a square lattice of size L 2 for L = 500. On 
Fig.l we plot the size Sk of the wave as a function of its 
number k in the avalanche. We can see that many of 
the next waves have the size greater than the size of the 
previous one. Moreover, even those waves, whose size is 
actually less than the size of its predecessor, are not most 
frequently embedded into the set of sites formed by the 
previous wave. For the particular avalanche presented on 
Fig. 1 the event of overlapping the previous wave by the 
next one occurs for all waves except the sixth and the 
last one. 

Anyway, one can note that a typical avalanche contains 
several sharp peaks corresponding to fast expansion of 
the avalanche size and in-between intervals of relatively 
slow, although irregular, contraction. One may expect 
that the average difference between subsequent waves in 
the slow phase follows a scaling law similar to Eq. ([!]). 
To verify this, one has to extract from the averaging of 
A(sfe) those waves which are related to the expansion 
phase. We can avoid this cumbersome and ambiguous 
procedure introducing new variables characterizing the 
'local' contraction and expansion. 

Consider two typical subsequent waves of topplings Wk 
and Wk+i with the sizes Sk and Sk+i, the (k + l)st wave 
overlaps the A:th wave. Let W be their intersection hav- 
ing the size s. 

Define the variables A + (sfe) = Sk — s and A~(sk) — 
Sfc+i — s, the first quantity is 'local contraction', the sec- 
ond one refers to 'local expansion'. We calculated the 
averages (A(s/ C )) for clusters, (A + (sk)) and (A~(sk)) 
for waves of topplings using data of 10 6 avalanches for 
the system size L = 500 (Fig. 2). The simulations show 
a power-law behavior (A(sfc)) ~ s a for clusters and 
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(A + (sfe)} ~ s a for waves; the exponents a and a + have 
close values. The value of the exponent a~ , for the rela- 
tion (A _ (sfc)) ~ s a is much smaller than a + . 

Concerning the estimation of the exponents a, a + and 
oT , we have to point out that the numerical determi- 
nation of these values is a rather difficult problem be- 
cause of a slow convergence of data obtained for large 
lattices to their limiting values. So, our numerical re- 
sults a w a+ w 0.88 and a~ w 0.29 for L = 500 are still 
far from the expected limit. The problem of estimation 
of these exponents is somehow similar to the numerical 
determination of the exponent r (for discussion, see, for 
example, [jL3[ ). The extrapolation L — > oo gives us some 
wide interval of possible values of a and a + which in- 
cludes 3/4, the theoretical prediction [1(J for the exact 
value. 

Being equivalent to r from a computational point of 
view, the exponent a + is more convenient for theoret- 
ical evaluations. The spanning tree representation of 
waves H makes it possible to interpret A + (s) as a sum 
of branches attached to the boundary of a wave 10 1, and 
then to use exact results obtained for the one-component 
Potts model @. 

The relative magnitude of A _ (sa,) and A + (sa,) is such 
that for large s the contraction of avalanche dominates 
its expansion. We show on Fig.l by filled diamonds the 
sizes of waves Sk subtracted by the size A~(sk-i)- It 
is clear that neglecting the quantities A~(sk) we do not 
change the qualitative dynamical picture of the avalanche 
and the contraction of waves can be described in terms 
of A + (sfc). Based on these data we can also see that the 
average (A(s w )} for waves which is equal to the remain- 
der (A+(sfc)) — (A~(sfc)) is actually negative for small 
waves and positive for large ones, as it was found in Jl2| . 

Finally, we establish the relationship of the type of 
Eq.(7) for the exponents r and a + . Following argumen- 
tation for clusters (section II), we estimate the known 
asymptotics of size distribution P(s w ) of waves Eq.(||) 
which is proportional to the probability of avalanches 

and to 



whose size S exceeds s„ 
the density of waves 

P{s w ) - 



P(S > s w ) ~ s v 



-T + l 



dn 



(8) 



Let N(s, s — t) be the number of waves in a particular 
avalanche with sizes between s — t and s provided that 
the size of the given avalanche is greater than s. The 
asymptotic behavior of dn/ds w can be evaluated as 



for large s v 
t where c = 
obtain 



dn 
ds w 

> t. 
0(1) 



(N(s w ,s w - 1)) 



t 



(9) 



Take in Eq.(g) c(A + {s w )) instead of 
is a constant. Then, from Eq.(0) we 



<N(s c ,s c -<A(s»> 

i <N(s w ,S w -2<A(s w )>)> 



FIG. 3. The expected number {N(s,s — t)) of clusters 
(waves) of size from the interval [s — t, s] provided that the 
size of the corresponding avalanche is greater than s. Due to 
fluctuations the data for waves are uncertain for some large 
values of s w , so the averages are not proportional. 

On Fig. 3 we present the quantity 

(N(s w ,s w -2(A+(s w )))) 

as a function of s w calculated from our numerical data 
where we put c = 2. It is apparently evident that asymp- 
totically it scales as O(l). Thus, we obtain from Eq.([lO|) 

T + a+ = 2. (11) 

For comparison on Fig. 3 we present the function 

(N(s c ,s c -(A(s c )})) 

for clusters which also scales as O(l) confirming Eq.(7). 



IV. ANALYSIS OF CONDITIONAL WAVE 
DISTRIBUTION 

Recently, Paczuski and Boettcher |l2j have undertaken 
careful numerical simulations to find the size distribution 
of subsequent waves for a given size of the preceding wave 
P(sk+i\sk)- The data was represented by a scaling form 



P( Sfc+1 | Sfe )~ Sfc -£F(^i), 



•Sfc 



(12) 



where F{x) — > const when x — > and F(x 3> 1) ~ x~ r . 
The function P(sk+\\sk) being considered as a normal- 
ized distribution of Sk+i should be multiplied by the fac- 
tor sf _1 to provide the normalization condition 



PN(s k +l\Sk)dSk^ 



x~ F(x)dx, (13) 



c(A + (s w )))). 



(10) 



where s co ~ L 2 is the cutoff in the wave sizes from 
the finite system size 12 1. The normalized function 
Pjv(sfe+i|sfe) has the asymptotics 
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TJ f I \ -0 0-1 

P N (s k+1 \S k ) ~ s fc+1 s^ , 



when s fe+ i <C Sfc, and 



-P/v(sfc+i|s fc ) 



-/3-r 0+r-l 
6 fe+l *fc i 



(14) 



(15) 



when s k+1 > s k . 

To be consistent with the analysis in Jl(]], both expo- 
nents /3 and r should be explained from the same point of 
view based on the spanning tree representation of waves. 
Before doing that, we will discuss an attempt to verify 
the existence of the scaling law Eq. ^ for the contraction 
phase of an avalanche by calculating the average differ- 
ence between subsequent waves (A(sfc)} = (s k — Sfc+i) 
p2[ . By the assumption |Tc| |, the average (A(sfc)} ~ 
s k >0 for the main part of an avalanche corresponding 
to the process of slow contraction of wave fronts. A seri- 
ous problem, however, is how to select the waves relating 
to the slow phase. In a real avalanche, at least waves cor- 
responding to the largest contribution to the expansion 
should be removed from the averaging as it was explained 
in Section II. Without the selection of waves responsible 
for the contraction of avalanches the result of averaging 
of As obtained in jl^] is easily predictable. 

Following Jl2| , fix a value s and take all waves with 
Sfe = s together with the subsequent waves of size s k +i 
from all avalanches whose size S > s. Consider sepa- 
rately the cases s^+i < Sk and s k +i > s k . In the first 
case, the average difference (As' - '} = (s k — s k +i) is ob- 
viously (As( - ') < s. 

In the opposite case, the waves with s k +i > s k have a 
power law asymptotics P(s k +i) ~ s k , 1 where 1 < 9 < 2 
for all s k . Indeed, the size distribution of waves with 
an arbitrary origin is P(s) ~ s^ 1 according to the two- 
component spanning tree representation The distri- 
bution of waves of size s with the origin in a fixed unique 
site is -P'(s) ~ s~ 2 . In the considered case, the subse- 
quent wave starts at a site in the localized area inside 
the previous wave. This implies 1 < 9 < 2. Therefore, 
the averaged positive difference {As^ + ') = {s k — s k+ i) ~ 
£2(2-6) _ s anc j diverges with the lattice size L. Thus, 
(As( + )) > (As(~)} until s becomes large and finite-size 
effects become essential. We see that the negative values 
of (As(sfc)) obtained in Q are not actually surprising 
and indicate only that the simple average (As) does not 
exhibit a power-law dependence on s k and cannot be re- 
lated directly to the density of waves. 

Nevertheless, the distribution Eq. (|l2|) itself brings im- 
portant information on avalanches. The exponents char- 
acterizing its asymptotics are related to basic exponents 
of the sandpile model. 



Consider first the exponent (3. According to Eq.(12), 
this exponent determines the behavior of smaller waves 
following just after waves of larger sizes: s k +\ <C s k . 
The wave W k can be represented by a tree T k covering 
the area s k and having the root i at the point where the 
wave W k was initiated ||. The tree T k+ \ representing 
the wave W k +i has the root at the same point i. 




FIG. 4. The structure of spanning trees representing two 
subsequent waves with sizes s k +i < s k - The origin of the 
avalanche is denoted by i. B is the point of intersection of 
the boundaries of these waves. 

To provide the sharp reduction of the next wave, T k 
and Tfe+i must have a special structure. Since each site in 
a wave topples exactly once, the state of the system inside 
a wave does not change after the wave is completed. This 
means that any branch of T k+ \ attached to the root i re- 
peats a branch of T k until it ends inside s k or touches the 
boundary of s k . It follows from s k+ i s k that at least 
one point B exists where the boundaries of W k and W k +\ 
intersect (Fig. 4). In the vicinity of B, both trees T k and 
Tfc+i can be attached by a bond b to their complemented 
subtrees T k and T k+1 || defined by the condition that a 

unification T and T gives a complete spanning tree of 
the whole lattice. Inversely, deletion of the bond b from 
the corresponding trees produces subtrees T k and T k+ i 
which obey the statistics of disconnected branches of lat- 
tice spanning trees [p~5|] or, equivalently, the asymptotics 
of last waves M 



Plast(s) 



1 



(16) 



It has been demonstrated in [Jl2[ that the function F(x) 
in Eq.([l2]) is constant in a finite interval < x < c < 1. 
Thus, we can consider the distribution 



P c (s k+ i) = J P N {s k+ i\s k )ds k 



(17) 



Sk>Sk + l/c 



with the function Pisf{sk+i\sk) taken in the form Eq.(pT[) 
This gives 



P c (s k+1 ) = s k ^ J sf 1 ds k 



s k >s k + 1 /c 



L 2P 

IT' (18) 
s k+i 
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On the other hand, P c (sk+i) apart from the normaliza- subtrees and Tk+i are strongly connected (Tk+i is a 
tion factor is given by a joined probability distribution of part of T^), the distributions of their sizes can roughly 
disconnected branches Tk and Tk+i- Despite the fact that be considered as independent. Then, we obtain 



Plast(sk+l)Plast(sk > Sfe+l/c) 



-11/ 
S fe+1 



-11/8 



dsi 



s k >s k + 1 /c 



7/4 ' 
S fe+1 



(19) 



To get P c (sk+i) from Eq.(|n]), we have to multiply the 
last expression by s^+i, the number of possible positions 
of the root i inside the disconnected branch Tk+i- Fi- 
nally, P c (sk+i) ~ V s fc+i an( i comparing with Eq. (|l8l ) 
we get (3 = 3/4: which explains the numerical result |12|j . 

To relate the exponent r with the exponent r in the 
distribution of a number of distinct sites covered by an 
avalanche, we consider waves of two types. A wave will be 
referred to as the growing one or the G-wave if Sj+i > 
and the reducing one or the R-wave if Sj+i < Sj. Every 
avalanche corresponds to a unique sequence of G-waves 
and R- waves e.g. GRGGR.... 

The number of distinct sites Sd in an avalanche is pro- 
portional to the size of the maximal wave W max , so we 
can expect that 



P\Smax) 



1 



(20) 



with the same critical exponent r. 

The expected number of waves in an avalanche diverges 
logarithmically with the lattice size ||. However, if the 
idea about the fast expansion phase is correct, the ex- 
pected number of waves in the interval between the maxi- 
mal wave and the latest wave Wk with Sk ~ 0(1) before 
the maximal wave should be finite when L — ► oo. 

Starting with this assumption, consider a finite se- 
quence of n G-waves and R-waves between the waves 
Wk +i and Wk + m = W max (for simplicity, we denote 
their sizes by Si, s n ). The first and the last waves in 
the sequence are clearly of type G. It follows from the 
numerical data of |L2] that the asymptotics Eqs.pl) and 



P(Sn) 



( |l5| ) of the distribution function P^jsk+i [sfc) are factor- 
ized. Extrapolating the distribution Eq.(fi~3) to the case 
of vanishing previous waves Sk — > 1, we can obtain the 
distribution of first waves in the sequence: 



1 



(21) 



For an avalanche GG... beginning from two G-waves, 
the distribution of the second wave is given by 



*2 

P(s 2 ) = J ds 1 P{s 1 )P N (s 2 \s 1 ). 



(22) 



Using Eqs.([l5|) and ( pl| ) we have for large s 2 

hlS2 (23) 



s 2 



Similarly, the leading asymptotics for the n-th wave in a 
sequence of n G-waves is 



P(s n ) 



(lns n y 



(24) 



The presence of R-waves reduces the logarithmic diver- 
gence of the numerator. For instance, using Eq. (|l4|) and 
Eq.(|l5|) we get in the case GRG... the numerator Insz 
instead of {Ins^j 2 in the case GGG... . Generally, if k 
(k < n — 2) last G-waves in the sequence follow R-wave, 
the convolution 



ds 1 ...ds n - 1 P(s 1 )P N (s 2 \s 1 )...P N (s 

n | Sn— 1 ) 



(25) 



has the asymptotics 



from scaling arguments. 



P{Sn) 



(lns n y 



(26) 



V. DISCUSSION 



Thus, for any finite sequence of G-waves and R-waves 
between the relatively small fco-th wave and the max- 
imal wave, we have the leading exponent [3 + r which 
governs the distribution of the maximal waves P max ~ 
s max ~ s max r - The numerical values obtained in are 
(3 = 3/4, r = 1/2. This gives r = 5/4 obtained in M] 



In conclusion, the analysis of the decomposition of 
avalanches into waves of topplings shows that a difference 
between two subsequent waves can be described by ap- 
propriate variables which follow a power law dependence 
on the wave size s. The exponent a + corresponding to 



G 



the contraction of waves can be related to one of the basic 
avalanche exponents r. 

The relation between the asymptotics Eq.(|l5|) of the 
distribution of subsequent waves -Pjv(sfc+i|sfc) in the case 
Sfe+i 3> s/c and the exponent t in the distribution of 
distinct sites involved into an avalanche implies an al- 
ternative way of determining r. Instead of derivation r 
from the analysis of slow contraction process, we can use 
the statistics of large waves Wk+i overlapping their pre- 
decessors Wk to link r with the exponents (3 and r. This 
approach sheds new light on the renormalization group 
(RG) procedure proposed by Pietronero et al. || for the 
sandpile model. In the RG method, one deals with sites 
of three classes: stable, critical and unstable. Extend- 
ing the characterization of the stationary properties at a 
generic scale, one describes the propagation of instabil- 
ity through the lattice taking into account only one-shot 
relaxation events at each scale. Thus, proliferation ef- 
fects due to multiple relaxations are not considered in 
this scheme. In this respect, the process described by 
RG is not a true avalanche, rather it is a wave propagat- 
ing from a given point or from a cluster of a given size. 
Correspondingly, the critical exponent determined in this 
way is actually the sum of exponents f3 + r in the asymp- 
totics of distribution of large waves. Its numerical value 
1.248 obtained in [|l6| is in excellent agreement with the 
value (3 + r — 5/4 proposed in |12|. On the other hand, 
it was shown in Section IV that (3 + r = r which explains 
the validity of the RG approach despite the neglect of 
multiple relaxations. 
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